Orbital fluctuations and strong correlations in quantum dots 



Gergely Zarand 
Theoretical Physics Department, Institute of Physics, 
Budapest Univ. of Technology and Economics, Budafoki ut 8, H-1521 
(Dated: February 2, 2008) 

In this lecture note we focus our attention to quantum dot systems where exotic strongly correlated 
behavior develops due to the presence of orbital or charge degrees of freedom. After giving a concise 
overview of the theory of transport and Kondo effect through a single electron transistor, we discuss 
how 5C/(4) Kondo effect develops in dots having orbitally degenerate states and in double dot 
systems, and then study the singlet-triplet transition in lateral quantum dots. Charge fluctuations 
and Matveev's mapping to the two-channel Kondo model in the vicinity of charge degeneracy point 
are also discussed. 
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I. INTRODUCTION 

Although no rigorous definition exists, we typically call 
a 'quantum dot' a small artificial structure containing 
conduction electrons, and weakly coupled to the rest of 
the world. There is a variety of ways to produce these 
structures: Maybe the most common technique to do 
this is by denning a typically /im-size region by shap- 
ing a two-dimensional electron gas using gate electrodes 
placed on the top of a semiconductor heterostructure or 
by etching (see, e.g., Refs. 1 2.3). In Fig. ^we show the 
top view of such a single electron transistor (SET) that 
has been first used to detect the Kondo effect in such 
a structured Beside semiconductor technologies, quan- 
tum dots can also be built from metallic grains^, and 
more recently it became possible to integrate even real 
molecules into electronic circuits^. The common feature 
of all these devices is that Coulomb correlations play an 
essential role in them, and induce Coulomb blockade- 
and Kondo effect iSui 

In the present paper (which has been prepared as a 




FIG. 1: Top view of the SET used by David Goldhaber- 
Gordon and his collaborators to first observe the Kondo ef- 
fect in a quantum dot, from Ref. Q: The white areas indicate 
regions where conduction electrons are present. The quan- 
tum dot is at the central region (white circle). The various 
electrodes (1-4) have been used to define the dot and the 
junctions. 



lecture note) we shall not attempt to give a complete 
overview of this enormous field. Instead, we shall first 
give a concise introduction into the basic properties of 
these devices, and then focus our attention to some ex- 
otic strongly correlated states associated with orbital and 
charge degrees of freedom that appear in them. 

There is two essential energy scales that characterize 
an isolated quantum dot: One of them is the charging en- 
ergy, Ec , the typical cost of putting an extra electron on 
the device. The other is the typical separation of single 
particle energies, also called level spacing, A. Typically 
A <C Ec, but for very small structures {e.g. in the ex- 
treme case of a molecule) these two energy scales can 
be of the same order of magnitude. While the charging 
energy Ec is usually of the order of e 2 /L with L the 
characteristic size of the device, the level spacing A de- 
pends very much on the material and dimensionality of 
the dot: it is typically very small in mesoscopic metallic 
grains, where it roughly scales as A ~ Ep/(kpL) 3 , and 
becomes of the order of one Kelvin only for nanoscale 
structures with L ~ 20A For two-dimensional semi- 
conductor structures, on the other hand, both kp and 
the Fermi energy are much smaller and since A scales as 
A ~ Ep J '(kpL) 2 , in these structures A becomes of the 
order of a Kelvin typically for L ~ 0.1/im. 

Coulomb correlations may become only important if 
the measurement temperature is less than the charging 
energy, T < Ec- Clearly, this criterion can be only sat- 
isfied with our current cooling technology if Ec is in the 
range of a few Kelvins, i.e. the size of the system is in 
the fj,m range or below. The behavior of a quantum dot 
is also very different in the regimes T > A and T < A: 
while in the former regime electron-hole excitations on 
the dot are important, for T < A these excitations do 
not play an essential role. 

Beside the difference in the typical energy scales Ec 
and A, there is also a difference in the way semiconduct- 
ing and metallic devices are usually connected. Although 
metallic particles can also be contacted through single or 
few mode contacts using e.g. STM tips, these grains 
are typically connected through multichannel leads with 
contact sizes much larger than the Fermi wavelength Xp. 
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Lateral semiconducting devices are, on the other hand, 
usually contacted through few or single mode contacts 
(though e.g. vertical dots are connected through a large 
contact area and thus many channels). While these de- 
tails can be important for some phenomena^S*ii the be- 
havior of all these devices is very similar in many re- 
spects. In the following, we shall therefore mainly focus 
on lateral quantum dots with single mode contacts. 

This lecture is organized as follows: First, in Sec- 
tion [n] we shall discuss the phenomenon of Coulomb 
blockade and the basic Hamiltonians that are used to de- 
scribe quantum dots. In Section IlIII we discuss how the 
Coulomb blockade is lifted by the formation of a strongly 
correlated Kondo state at very low temperatures, T <C A. 
In Section llVI we shall study more exotic strongly corre- 
lated states that appear due to orbital degrees of free- 
dom, the SU(4) Kondo state, and the so-called singlet- 
triplet transition. Section E] is devoted to the analysis 
of the Coulomb blockade staircase in the vicinity of the 
degeneracy points, where strong charge fluctuations are 
present. 



II. COULOMB BLOCKADE 

In almost all systems discussed in the introduction we 
can describe the isolated dot by the following second 
quantized Hamiltonian 



Hdot = ^2 e 3 
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where the second term describes interactions between the 
conduction electrons on the island, and the effects of vari- 
ous gate voltages are accounted for by the last term. The 
operator d^ a creates a conduction electron in a single par- 
ticle state tfj with spin a on the dot. 

Fortunately, the terms He = H- mi + H gatc can be re- 
placed in most cases to a very good accuracy by a simple 
classical interaction termiS 



Hc = 2C ' " dot 



v„c g 



(2) 



where C denotes the total capacitance of the dot, C g is 
the gate capacitance, e is the electron's charge, V g stands 
for the gate voltage (roughly proportional to the voltage 
on electrode '2' in Fig. 0), and Udot = Ylj a • dj a dj a is 
the number of extra electrons on the dot. This simple 
form can be derived by estimating the Coulomb integrals 
for a chaotic doti£, however, it also follows from the phe- 
nomenon of screening in a metallic particle. We outline 
this rather instructive derivation of Eq. (J2J within the 
Hartree approximation in Appendix ^ Clearly, the di- 
mensionless gate voltage N g — V g C g /e sets the number 
of electrons on the dot, (ridot) ~ N g . 

The single particle levels e,- above are random but cor- 
related: The distribution of these levels for typical (i.e. 
large and chaotic) islands is given with a good accuracy 
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FIG. 2: Excitation spectrum of an island. Lines represent 
eigenenergies of the island. Charging excitations typically 
need an energy ~ Ec while internal electron-hole and spin 
excitations cost an energy A <C Ec- 



by random matrix theory^ which predicts among oth- 
ers that the separation s between two neighboring states 
displays a universal distribution, 



P( s ) = ip( s /A) 



(3) 



with A the average level spacing between two neighbor- 
ing levels. For small separations, energy levels repel each- 
other and p vanishes asp(^) ~ (x)^) where the expo- 
nent is (3 — 1, 2 or 4, depending on the symmetry of the 
Hamiltonian (orthogonal, unitary, and simplectic, respec- 
tively) . In some special cases cross-overs between various 
universality classes can also occur, and in some cases level 
repulsion may be even absent for dots with special sym- 
metry properties. 

The spectrum of an isolated dot described by Eqs. (|TJ 
and J5J is sketched in Fig. As we already mentioned, 
for typical parameters and relatively large lateral dots 
or metallic islands, the charging energy is much larger 
than the level spacing, Ec S> A. Accordingly, internal 
electron-hole excitations cost much smaller energy than 
charge excitations of the dot. For dot sizes in the 0.1/J,m 
range the capacitance C can be small enough so that the 
charging energy Ec — e 2 /2C associated with putting an 
extra electron on the dot can safely be in the meV range. 



Therefore, unless N a 



VaC 



is a half-integer, it costs 



a finite energy to charge the device, and therefore the 
number of electrons on the dot becomes quantized at low 
enough temperatures, T -C Ec, and a Coulomb blockade 
develops - provided that quantum fluctuations induced 
by coupling the dot to leads are not very strong. 

Let us now consider a quantum dot that is weakly 
tunnel-coupled to leads, 'weak coupling' in this context 
meaning that the conductance between the island and the 
leads is less than the quantum conductance, Gq = 2e 2 /h. 
In the particular case of a lateral quantum dot this con- 
dition is satisfied when the last conduction electron chan- 
nel is being pinched off. If the lead-dot conductance is 
much larger than Gq then the Coulomb blockade is lifted 
by quantum fluctuations (see also Section 0) . In the 
absence of tunneling processes, the charge on the dot 




FIG. 3: The number of electrons on a quantum dot as a func- 
tion of the dimensionless gate voltage. The sudden jumps of 
an isolated dot become smeared out due to quantum fluctua- 
tion as soon as we couple the dot to leads. 



would change in sudden steps at T = temperature (see 
Fig. I2J. However, in the vicinity of the jumps, where 



V a C a 



is a half-integer, two charging states of the is- 
land become almost degenerate. Therefore quantum tun- 
neling to the leads induces quantum fluctuations between 
these two charging states, smears out the steps, and 
eventually completely suppresses the steps only leaving 
some small oscillations on the top of a linear (ndot)(^g) 
dependence^ At a finite temperature T / 0, thermal 
fluctuations play a similar role, and the charging steps 
are also washed out if T ^> Ec- 

Charge quantization is also reflected in the transport 
properties of the dot. To study transport through a quan- 
tum dot, one typically builds a single electron transistor 
(SET) by attaching the dot to two leads, as shown in 
Figs. ^ and 0] Let us first assume that the conductances 
Gl and Gr between the dot and the lead on the left / 
right are small compared to Gq and that quantum fluc- 
tuations are small. In this limit we can describe charge 
fluctuations on the dot by the following simple tunneling 
Hamiltonian: 
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(4) 



where we assumed single mode contacts. The fields i$ L e(T 

and tp R ecr denote the creation operators of conduction 
electrons of energy e and spin a in the left and right leads, 
respectively, and are normalized to satisfy the anticom- 
mutation relations {w L / R ea i ipL/R,e'cr'} = S acr 'S ee '. Note 
that the tunneling matrix elements tj and tj fluctuate 
from level to level since they depend on the amplitude of 
the wave function tpj at the tunneling position. Using a 
simple Boltzmann equation approach one then finds that 
for T < Ec the linear conductance of the SET has peaks 
whenever two charge states of the dot are degenerate and 
N g is a half-integer. In the regime, A <C T <C Ec, in par- 
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FIG. 4: (a) Sketch of a single electron transistor, (b) 
Conductance of the SET as a function of gate voltage for 
Ec <C T <C A. At the degeneracy points Coulomb blockade 
is lifted and transport is allowed through the single electron 
transistor. 



ticular, one finds thati^ 

G(T)«i GlGr 



AE/T 



2G L + G r smh(AE/T) 



(5) 



where AE is the energy difference between the 
two charging states of the dot, and Gl/r = 

(8ir 2 e 2 /h) g^ot Ql/r (\tf \ 2 )j is the tunnel conductance 
of the two junctions, with gdot and Ql/r the density of 
states on the dot and the leads. Note that even for per- 
fect charge degeneracy, AE = 0, the resistance of the 
SET is twice as large as the sum of the two junction resis- 
tances due to Coulomb correlations. As we decrease the 
temperature, the conductance peaks become sharper and 
sharper, while the conductance between the peaks de- 
creases exponentially, and the Coulomb blockade devel- 
ops. This simple Boltzmann equation picture, however, 
breaks down at somewhat lower temperatures, where 
higher order processes and quantum fluctuations become 
important. These quantum fluctuations may even com- 
pletely lift the Coulomb blockade and result in a perfect 
conductance at low temperatures as we shall explain in 
the next section. 

Let us first study the conductance of the SET in the 
regime where the difference AE between the energy of 
the two charging states considered is much larger than 
the temperature. It turns out that the range of valid- 
ity of Eq. JSJ describing an activated behavior is rather 
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Fig.a 
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Fig.b 

FIG. 5: Elastic (a) and inelastic (b) co-tunneling processes. 
Inelastic co-tunneling processes give a conductance ~ T 2 
while elastic co-tunneling gives a finite conductance as T — > 0. 



small for typical parameters, and the conductance is 
dominated by second order virtual processes as soon as 
we lower the temperature much below AE. From the 
point of view of these second order processes two regimes 
must be distinguished: In the regime A <C T <C Ec 
the leading term to the conductance comes from elastic 
and inelastic co-tunneling processes shown in Fig. [21 16 
In the inelastic co-tunneling process a conduction elec- 
tron jumps into the dot from one lead and another elec- 
tron jumps out of the dot to the other lead in a sec- 
ond order virtual process, leaving behind (or absorb- 
ing) an electron- hole excitation on the dot, while in an 
elastic co-tunneling it is the same electron that jumps 
out. Inelastic co-tunneling gives a conductance Gi nc i ~ 
GlGr{T / Eq) 2 /(e 2 /h), and is thus clearly suppressed as 
the temperature decreases^ while elastic co-tunneling 
results in a small temperature independent residual con- 
ductance, G c i ~ G L G R (A/E c )/(e 2 /h) even at T = 
temperature^ 



left right 




FIG. 6: Exchange process leading to the enhancement of the 
conductance as T — > 0. 



III. KONDO EFFECT 

For T -C A <C Ec inelastic co-tunneling processes are 
not allowed, and the properties of the SET depend essen- 
tially on the number of electrons on the dot. The ground 
state of the isolated dot must be spin degenerate if there 
is an odd number of electrons on the dot, while it is usu- 
ally non-degenerate, if the number of electrons on the 
dot is even. In the latter case nothing special happens: 
quantum fluctuations due to coupling to the leads pro- 
duce just a residual conductance as T — > 0. If, however, 
the number of the electrons on the dot is odd, then the 
ground state has a spin degeneracy, which can give rise to 
the Kondo effect discussed below. In this case, exchange 
processes shown in Fig. [5] give a contribution to the con- 
ductance. As we lower the temperature, the effective 
amplitude of these processes increases due to the Kondo 
effect, and ultimately gives a conductance that can be as 
large as the quantum conductance Gq = 2e 2 /h at T = 
temperature^ This strong enhancement is due to strong 
quantum fluctuations of the spin of the dot, and the for- 
mation of a strongly correlated Kondo state. The typical 
temperature dependence of the conductance for N g odd 
is shown in Fig. [7] 
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FIG. 7: (a) Linear conductance of the SET for T K < T < 
Ec (dashed line) and T — > (solid line), (b) Temperature 
dependence of the conductance of a symmetrical SET for an 
odd number of electrons on the dot (continuous line) and for 
an even number of electrons with non-degenerate ground state 
(dashed line). 

To understand why the conductance of the dot be- 
comes large, let us keep only the last, singly occupied 
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level that gives rise to the Kondo effect, dj a 
€j — ► £rf, and write the Hamiltonian of the dot as 

H = J2 e ddUa + H c 

a 

+ E E / dee^l a (e)^ >a (e) 

~ T D J 



(6) 



a p 1 =L,R 

de \t L g 1 ^ 2 (4^,a(e) + h.c.)+ " 



1 /2 

where we introduced the fields, ipL,a{^) = q//r ^L/R.ea- 
For the sake of simplicity, let us assume that ^g 1 ^ 2 = 

T 1/2 

t g R ■ If we then make a unitary transformation and 
introduce the even and odd field operators, 



V2 



(7) 



then obviously the odd combination ip_ fully decouples 
from the dot, and the tunneling part of the Hamiltonian 
can be written as 



V = t de (4v+ )CT (e)+h.c.) 



(8) 



where t — \/2 t R g R ^ ■ One can now perform second or- 
der degenerate perturbation theory in V in the subspace 
J2 a d]jd a = 1 to obtain the following effective exchange 
Hamiltonian: 



H, 



e^+AM+Ae) ( 9 ) 
de de'^ +a {e)a aal ip +IJ ,{e') , 
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where S — ^(d)ad) is the spin of the dot, and J ~ t 2 /Ec 
is a dimensionless antiferromagnetic exchange coupling. 
Thus electrons in the even channel ip + couple antiferro- 
magnetically to the spin on the partially occupied c£-level, 
and try to screen it to get rid of the residual entropy as- 
sociated with it. 

Written in the original basis, Eq. Q contains terms 
~ J Sip R (e)aipL(i') , which allow for charge transfer from 
one side of the dot to the other side, and in leading order 
these terms give a conductance ~ J 2 . However, higher 
order terms in J turn out to give logarithmically diver- 
gent contributions, 



G (j 2 + 2 J 3 ln(A/T) 
h V 



(10) 



As a result, the conductance of the device increases as we 
decrease the temperature and our perturbative approach 
breaks down at the so-called Kondo temperature, 



-1/j 



(11) 



One can try to get rid of the logarithmic singularity in 
Eq. l|lUfl by summing up the most singular contributions 



in each order in J. This can be most easily done by 
performing a renormalization group calculation and re- 
placing J by its renormalized value in the perturbative 
expression, G ~ ^-J 2 .— However, this procedure does 
not cure the problem and gives a conductance that still 
diverges at T = Tk, 



G 



1 



h \n 2 (T/T K ) 



(12) 



The meaning of the energy scale Tk is that below this 
temperature scale the effective exchange coupling be- 
comes large and a conduction electron spin is tied anti- 
ferromagnetically to the spin of the dot to form a singlet 
(see Fig.®. 



It is not difficult to show that then for t L g 1 / 2 



the SET must have perfect conductance at T = temper- 
ature. To show this, let us apply the Friedel sum rulejiiiS 
that relates the number iVbound of electrons bound to the 
impurity and the phase shifts S of the electrons ip+ 1 as 



Abound — 2— , 
7T 



(13) 



where the factor 2 is due to the spin. This relation 
implies that in the even channel conduction electrons 
acquire a phase shift S = tt/2, and correspondingly, 



04 



, 2 it 



04 



-ip + in course of a scattering process at 



T = temperature. Going back to the original left-right 
basis, this implies that left and right electrons scatter as 



-i>i 



(14) 



In other words, an electron coming from the left is trans- 
mitted to the right without any backs cattering, and thus 
the quantum dot has a perfect conductance, 2e 2 /h. 

We remark here that only a symmetrical device can 
have a perfect transmission, and only if the the number 
of electrons on the d-level is approximately one, (rid) ~ 1- 
All the considerations above can be easily generalized to 



the case t L g^ ^ t R 6R i an d one obtains for the zero 
temperature conductance in the Kondo limit, 
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(15) 



which is clearly less than 2e 2 /h for non-symmetrical dots. 

The phase shift 6 — tt/2 also implies that there must be 
a resonance at the Fermi energy. In fact, this resonance 
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FIG. 8: Conduction electrons in the even channel screen the 
impurity spin. 
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FIG. 9: The Kondo resonance appears as a peak in the differ- 
ential conductance of the SET (figure taken from Ref. 0). At 
T = temperature the conductance approaches the quantum 
limit of the conductance, 2e 2 /h. 

is called the Kondo resonance, and can be directly seen 
in the differential conductance (related to the density of 
states as usual) of the single electron transistor shown in 
Fig. EjJ This is a many-body resonance that develops at 
the Fermi energy (zero bias) as the temperature is cooled 
down below the Kondo temperature Tk- 

The basic transport properties of the SET have been 
summarized in Fig. [7] Although we could get a fairly 
good analytical understanding of behavior of a SET, 
based on the simple considerations outlined above, it is 
rather difficult to obtain a quantitative description. In 
fact, to obtain a quantitative description extensive com- 
putations such as numerical renormalization group cal- 
culations are needed 2 ^. A valid and complete descrip- 
tion of the out of equilibrium physics of a SET is still 
missing i 21 i 22 i 23 i 24 

IV. ORBITAL DEGENERACY AND 
CORRELATIONS IN QUANTUM DOTS 

In the previous section we sketched the generic behav- 
ior of a quantum dot, and assumed that the separation of 
the last, partially occupied level is by an energy distance 
~ A 3> T separated from all other single particle levels 
of the dot. This is, however, not always true. In the 
case of a symmetrical arrangement like the one shown 
in Fig. El e.g., some of the states of the dot are or- 
bitally degenerate by symmetryr^L while in some other 
cases almost degenerate orbitals may show up just by 
accident i^SiSS This orbital degeneracy can play a very im- 
portant role when we fill up these degenerate (or almost 
degenerate) levels, and leads to such phenomena as the 
SU(4) Kondo stat e 25 ' 27 or eventually the singlet-triplet 
transitioni 26 i 28 i 29 

SU(4) physics and triangular dots 

Let us first discuss what happens if we have just a sin- 
gle electron on an orbitally degenerate level. The pres- 
ence of the orbital degeneracy leads to an unusual state of 
(approximate) SU(4) symmetry in this case, where spin 
and charge degrees of freedom are entangled. This SU(4) 



FIG. 10: Arrangement with triangular symmetry and the 
structure of the four-fold degenerate ground state of the iso- 
lated triangular dot. 



state has been proposed first theoretically to appear in 
double dot systems and quantum dots with triangular 
symmetry in Refs. |53 and l2Fjl However, while there 
is no unambiguous experimental evidence of an SU(4) 
Kondo effect in double dot devices^ the SU(4) state 
has been recently observed in vertical dots of cylindrical 
symmetry^! as well as in carbon nanotube single electron 
transistors^ In both chases the degeneracy index is due 
to a chiral symmetry just as in Ref.l25i Other realizations 
of states with SU(4) symmetry have been also proposed 
later in more complicated systemst^iiM and also in the 
context of heavy fermions. 35 

For the sake of concreteness, we shall focus here to the 
case of the triangular dot shown in Fig.^J However, our 
discussions carry over with trivial modifications to the 
previously mentioned experimental systems in Refs. 31 
and l32t Let us first assume that we have two orbitally 
degenerate levels |±) that can be labeled by some chiral- 
ity index r = ±, and let us focus on the charging of this 
multiplet only. At the Hartree-Fock level, these levels of 
the isolated dot can be described by: 

ffdot - 4AEtt> + AE 8 TT ,)d T , a ~ ^-S 2 

r,r' ,<j 

+ ^ K+n _)2 + ^ (n+ _ n _ ) 2 ) (16) 

where d\ a creates an electron on the dot within the 
degenerate multiplet with spin a and orbital label r. 
The energy shift AE above is proportional to the 
(symmetrically-applied) gate voltage and controls the 
charge on the dot, while E TT i accounts for the split- 
ting generated by deviations from perfect triangular sym- 
metry (^2 T E TT — 0). We denote the total number 
of electrons in state r = ± by n T = d\ a d Ta: and 
S = ^ X)t a a' dt,j<Jcrcr'dra' is the total spin of the dot. 
The terms proportional to Eq and Eq are generated by 
the Hartree interaction, while that proportional to Jh in 
Eq. ((To)) is the Hund's rule coupling, generated by ex- 
change. This term has no importance if there is only a 
single electron on the dot. 

Let us first consider the linear conductance. The cur- 
rents Ij between leads j and the dot are related to the 
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FIG. 11: Structure of the Coulomb blockade peaks of the tri- 
angular dot. The multiplets are labeled by the corresponding 
irreducible representations. The arrows indicate the direction 
the peaks move when applying a Zeeman field. 



voltages Vj applied on them by the conductance tensor, 



j' 



(17) 



which further simplifies to Gjj> — ^G6jj>—G/2 for a sym- 
metrical system. A schematic plot of the conductance G 
as a function of AE is shown in Fig. ^2 The arrange- 
ment of the four Coulomb blockade peaks associated to 
the fourfold degenerate state is symmetrical, and the 
height of the four peaks turns out to be numerically al- 
most identical at high temperatures^ 

The most interesting regime in Fig. llll appears between 
the first two peaks. Here there is one electron on the dot, 
and correspondingly the ground state of the isolated dot 
is fourfold degenerate. Let us now tunnel couple the dots 
to conduction electrons in the three leads, ipj (j = 1, 2, 3). 
To simplify the Hamiltonian we first introduce new fields 
that transform with the same symmetry as the states |r), 



i>j -> fa 



1 

7! 



±i 27TJ/3 



(18) 



With this notation the tunneling Hamiltonian of a per- 
fectly symmetrical dot becomes 



V" 



tJ2 Jde (4„lMe) 



h.c. 



(19) 



To describe the fourfold degenerate ground state of 
the dot, we can introduce the spin and and orbital spin 
operators S and f , with S z = ±1/2 and T z = ±1/2 cor- 
responding to the states a = ± and r = ±. Likewise, 
we can introduce spin and orbital spin operators a and 
t also for the conduction electrons, and then proceed as 
in the previous section to generate an effective Hamil- 
tonian by performing second order perturbation theory 
in the tunneling V, Eq. l|19fl . The resulting interaction 
Hamiltonian is rather complex and contains all kinds of 
orbital and spin couplings of the type ~ T + r~Sa or 





FIG. 12: Example of a virtual process generating a coupling 
between the spin and the orbital spin. 



~ T z t z m-i— These terms are generated by processes like 
the one shown in Fig. 1121 and clearly couple spin and 
orbital fluctuations to each-other. 

Fortunately, a renormalization group analysis reveals 
that at low temperatures the various couplings become 
equal, and the Hamiltonian can be simply replaced by the 
following remarkably simple SU(4) symmetrical effective 
Hamiltonian (Coqblin-Schrieffer model), 



H eS (T -» 0) 



a, /3=1,.. 



€^fs 1/3) H 



(20) 



where the index a labels the four combinations of pos- 
sible spin and pseudospin indices, and the |a)'s denote 
the four states of the dot. The dynamical generation of 
this SU(4) symmetry can also be verified by solving the 
original complicated Hamiltonian by the powerful ma- 
chinery of numerical renormalization groupS*^ We re- 
mark here that the structure of the fixed point Hamilto- 
nian, Eq. (121 )[l is rather robust. Even if the system does 
not have a perfect triangular (or chiral) symmetry, the 
exchange part of the effective Hamiltonian at low tem- 
peratures will take the form Eq. I|20|) . and the effect of 
imperfect symmetry only generates some splitting E TT > 
for the orbitally degenerate levels and some small poten- 
tial scattering^ These terms, of course, break the SU(4) 
symmetry of Eq. H20|) . but represent only marginal per- 
turbations, and do not influence the physical properties 
of the system in an esential way if they are small. In a 
similar way, the SU(4) symmetric fixed point discussed 
here may be relevant even for systems with (approxi- 
mate) accidental degeneracy even if they do not have a 
perfect SU(4) symmetry. 

Similar spin and orbital entangled states apparently 
also show up in molecular clusters, but there they may 
lead to the appearance of unusual non-Fermi liquid 
statesi^S^i 

The Hamiltonian Eq. (|20() is one of the exactly solvable 
models^ and has been studied thoroughly beforei^Llil 
Just as in the Kondo problem studied in the previous 
section, the SU(4) spin of the dot is screened below the 
'SU(4)' Kondo temperature, T^' . However, to screen 
an SU(4) spin, one needs three conduction electrons, as 
shown in Fig. As a result, the Friedel sum rule in the 
present case is modified to 



3 = ^^1=4^ 



7T 



7T 



(21) 
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impurity electrons 1 — 1 

FIG. 13: To screen the SU(4) spin of a triangular dot one 
needs three conduction electrons. The singlet formed corre- 
sponds to the Young tableau on the right, while the defin- 
ing four- dimensional SU(4) representations are denoted by 
squares. 



corresponding to a phase shift 6 = Sir/ '4 for the electrons 

The application of a magnetic field on the dot, H — > 
H — BS Z clearly suppresses spin fluctuations!* 1 - However, 
it does not suppress orbital fluctuations, which still lead 
to a more conventional SU(2) Kondo effect, by replac- 
ing the spin in the original Kondo problem. The Kondo 
temperature of this orbital Kondo effect is, how- 

ever, somewhat reduced compared to the SU(4) Kondo 



temperature j& 



T, 



K 



(T<?>) 



(22) 



The phase shifts in this case are simply <5j = 5±^ rs 
7r/2, while Si = 6± ± ~ 0, just as for the original Kondo 
problem. The splitting of the two levels r = ± has a 
similar effect and drives the dot to a simple spin SU(2) 
state. 

The zero-temperature phase shifts can be related to 
the transport properties of the device: From the T = 
phase shifts S Ta one can construct the conduction elec- 
trons' scattering matrix in the original basis ipj and 
compute all transport coefficients using the Landauer- 
Buttiker formula. 25 i 42 i 43 The T = conductance G turns 
out to be independent of the magnetic field, and both for 
the SU(4) and orbital Kondo states one finds the same 
value, 

G( s) = ^£- 3 MB)) = g. 



9h 

a 

The 'polarization of the current, 

sin 2 ((5 T ) - sin 2 (^) 



P = 



sin 2 ((5 T ) -rsin 2 ^) 



(23) 



however, does depend on the magnetic field and takes the 
values P — and P = 1 in the SU(4) and orbital Kondo 
states, respectively. The phase shifts in Eq. (|23fl can be 
extracted with very high accuracy from the finite size 
spectrum computed via the numerical renormalization 
group procedure ; 27 ^ 6 ' 44 and the results (originally com- 
puted for the double dot system in Ref . l27j) are shown in 
Fig. [21 Clearly, this device can be used as a spin fil- 
ter: Applying a Zeeman field one can induce a large spin 



polarization at low temperatures while having a large 
~ e 2 jh conductance through the device. A slightly mod- 
ified version of this spin filter has indeed been realized 
in Ref. |32|, where two orbital states originating from dif- 
ferent multiplets have been used to generate the orbital 
Kondo i 




FIG. 14: Spin polarization of the current through the trian- 
gular dot as a function of magnetic field. 

To close the analysis of the SU(4) Kondo effect, let us 
shortly discuss how the SU(4) state emerges in capaci- 
tively coupled quantum dots, where it actually has been 
identified first. In this device, shown in Fig. 1151 the ca- 
pacitively coupled dots can be described by the following 
simple Hamiltonian 



H, 



dot 



^±( n+ -N g+ Y + ^-{n_-N g _Y+E c n + n_ 



where the dimensionless gate voltages N g ± set the num- 
ber n + and n_ of the electrons on the left and right dots, 
respectively. The last term is due to the capacitive cou- 
pling between the two dots, and it is essentially this term 
which is responsible for the SU(4) physics discussed. As 
shown in Fig. 1151 b. in the parameter space of the two-dot 
regions appear, where the two states (n+,n_) = (1,0) 
and (n + ,nJ) = (0,1) are almost degenerate, while the 
states (ri + ,n_) = (0,0) and (n+,n_) = (1,1) are pushed 
to higher energies of order ~ Ec- In the simplest, how- 
ever, most frequent case the states (1, 0) and (0, 1) have 
both spin S — 1/2, associated with the extra electron on 
the dots. Therefore, in the regime above for temperatures 
below the charging energy Ec and the level spacing A 
of the dots, the dynamics of the double dot is essentially 
restricted to the subspace {S z = ±1/2; n + — n_ = ±1}, 
and we can describe its charge fluctuations in terms of 
the orbital pseudospin T z = (n + — n-)/2 = ±5. Cou- 
pling the two dot system to leads, we arrive at the very 
same Hamiltonian as for the triangular dots, although 
with very different parameters. Much of the previous 
discussions apply to this system as well which, in ad- 
dition to being a good spin-filter, also exhibits a giant 
magneto- resistance 
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Singlet-triplet transition 
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(a) 



(0,1) 



(i,i) 



(0,0) 



1/2 



(1,0) 



(b) 



3/2 



(2,0) 



FIG. 15: (a) A virtual process leading to entanglement be- 
tween charge and spin fluctuations and the SU(4) Kondo state 
in the double dot device, (b) Charging states of the double 
dot device as a function of the dimensionless gate voltages 
N g ±. The colored region indicates the regime where the two 
states (1,0) and (0, 1) are almost degenerate. 



So far we discussed the case where there is a single 
conduction electron on the degenerate levels. The regime 
between the two middle peaks in Fig. 1111 where there is 
two electrons on the (almost) degenerate multiplet is, 
however, also extremely interesting. 

In this regime the Hund's rule coupling Jh in Eq. H16|) 
is very important. This coupling is typically smaller than 
the level spacing A in usual quantum dots. If, however, 
it is larger than the separation between the last occupied 
and first empty levels, <5e, then it gives rise to a triplet 
ground state with 5=1. In particular, such a spin S = 1 
state forms if one puts two electrons on the orbitally de- 
generate level of a symmetrical quantum dot discussed 
beforet^iS but almost degenerate states may also oc- 
cur in usual single electron transistors just by accident^ 
Since in many cases one can shift the levels and thus tune 
5e by an external magnetic fieldSSiS or simply by chang- 
ing the shape of the dot by gate electrodes^ one can 
actually drive a quantum dot from a triplet to a singlet 
state as illustrated in Fig. 



8e «J, 



t I 

8e » J , 



8e 



FIG. 16: The state of a quantum dot changes from a triplet 
to a singlet as the separation 5e between the last occupied 
and first empty level increases. 



While this transition has been first studied in vertical 
dotsjS^£^£i here we shall focus on the usual lateral ar- 
rangement of Fig. 1161 which has very different transport 



properties 



20.44.49.50 



For the sake of simplicity let us as- 



sume that we have a completely symmetrical device and 
that the two levels |±) participating in the formation of 
the triplet state are even and odd. In this case we can in- 
troduce the even and odd fields tp± = (i/)r ± iPl)/V% by 
Eq. J7J , which by symmetry only couple to the even and 
odd states, |+) and |— ), respectively. The hybridization 
term in this case simply reads 

V = t+Y, [de U\^ +cr {e)+hA 

Y J jde{dl rr ^_ <7 {e) + h.c). (25) 



+ t. 



To describe the isolated dot we can use the following 
simplified version of Eq. 11611 , 



■yS +— (n++n_) 



(26) 



It is instructive to study the triplet state of the dot 
first, 6e <C Jh- Second order perturbation theory in 
Eq. d2HI) in this regime gives the following Hamiltonian 
replacing © 



#<*=£££ e4 ff (e)W, CT (e) 

r— ± a e 

J+ 



(27) 



Sj2j j dede'^ +a {t)o aa ^ +<7 ,{e') 

cr.tx' 

S £ / / de rfe Vl CT (e)o ! (r<r'V'-<7'(e / ) 



Clearly, the even and odd electrons couple with different 
exchange couplings J± to the spin. However, now S is 
a spin S = 1 operator, and to screen it completely one 
needs to bind two conduction electrons to it. This implies 
that an electron from both the even and the odd chan- 
nels will be bound to the spin, and correspondingly two 
consecutive Kondo effects will take place at temperatures 



Ae 



> T_ w A e 



(28) 



This also implies that the conductance at T — tempera- 
ture must vanish in the Kondo limit by the following sim- 
ple argument^ Again, we can use the Friedel sum rule 
to obtain the T = temperature phase shifts 5± = it/2 
in both the even and odd channels. In the original left- 
right basis this implies that the lead electrons scatter 
as "0L l r ~ * ~i ) L/R^ *- e -i their wave function vanishes at 
the dot position (by Pauli principle), and ipL/R are com- 
pletely reflected. 
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It is a simple matter to express the T = temperature 
conductance in terms of the T — temperature phase 
shifts by means of the Landauer-Buttiker formula as^ 



e 2 ^ 



r,' = ^^sm 2 (5 +)<T -<5_ )CT ) 



(29) 



This formula immediately implies that the conduc- 
tance as a function of a Zeeman field B must be non- 
monotonic^ As we argued before, the conductance of 
the dot is small for B = at T = temperature. 
For T_ <C B <C T+, however, the Kondo effect in the 
odd channel is suppressed and correspondingly the phase 
shifts in this channel are approximately given by <5_f ~ tt 
and 5- i ~ 0, while the phase shifts in the even channel 
are still 8- a ~ tt/2, and thus by Eq. (|29|) the conduc- 
tance must be close to 2e 2 /h. For even larger magnetic 
fields, B > T + ,T_, the Kondo effect is killed in both 
channels, and correspondingly <5±| ~ ir and 5±i w 0, 
resulting in a small conductance again. The magnetic 
field dependence of the phase shifts obtained from a nu- 
merical renormalization group calculation and the corre- 
sponding conductance are shown in Fig. 1171 44 By general 
arguments, 49 similar non-monotonic behavior must occur 
in the temperature- and bias-dependence of the conduc- 
tance, as it has indeed been observed experimentallyiSSiiS 




FIG. 17: Phase shifts (top) and dip structure in the conduc- 
tance as a function of Zeeman field B on the triplet side of 
the transition at zero temperature, compute d by the numeri- 
cal renormalizaion group method (from Ref. 44). 



it reappears on the singlet side. However, there the size 
of the dip is not determined by the smaller Kondo scale, 
T_, but rather by the excitation energy of the triplet, 
~ Se — Jh i 26 i 44 ' 50 Note that the transition between the 
triplet and singlet states is smooth and the singlet-triplet 
transition is rather a cross-over than a phase transition 
in the above scenario. 




FIG. 18: Phase shifts (top) and corresponding conductance 
(bottom) as a function of orbital splitting Se at temperature 
T = (from Ref.0). 

Also, the picture outlined above changes substantially 
if both states |±) happen to have the same parity, and 
couple only to one of the fields ip+. In this case the 
conductance is small on the singlet side of the tran- 
sition and exhibits a dip as a function of tempera- 
ture/magnetic field. 50 However, the spin of the dot can- 
not be screened on the triplet side even at T = tem- 
perature, and a real Kosterlitz-Thouless type quantum 
phase transition occurs, where the T = temperature 
conductance has a jump at the transition point^ - The 
triplet phase in this case is also anomalous jS - and in fact 
is of a 'marginal Fermi liquid type'^i^ since the spin 
is not fully screened^SS Correspondingly, the conduc- 
tance saturates very slowly, and behaves asymptotically 
as G ~ est — 1/ In 2 (T + /T). Similar behavior is expected 
to occur if the smaller Kondo temperature T_ is much be- 
low the measurement temperature, and indeed a behav- 
ior i n ag reement with the Kosterlitz-Thouless scenario of 
Ref. |5(J has been observed in some experiments^ 



Clearly, the T = temperature conductance must be 
also small on the singlet side of the transition, Se 3> Jh, 
where the dot is in a singlet state and no Kondo effect 
occurs. However, in the vicinity of the degeneracy point, 
5e w Jh, the triplet and the singlet states of the dot are 
almost degenerate, and quantum fluctuations between 
these four states generate another type of strongly corre- 
lated state with an increased Kondo temperature and a 
large conductancei^S^^^ 

To compute the full T = conductance as a func- 
tion of Se non-perturbative methods such as numerical 
renormalization group are needed^ In the vicinity of 
the transition point the conductance goes up to 2e 2 /h, 
in perfect agreement with the experimental observations. 
While the non-monotonic behavior characteristic of the 
triplet state disappears in the vicinity of the transition, 



V. CHARGE FLUCTUATIONS AND 
TWO-CHANNEL KONDO EFFECT AT THE 
DEGENERACY POINT 

In the previous section we focussed our attention to 
the regimes where charge fluctuations of the dot were 
negligible. In the vicinity of the degeneracy points, N g sa 
half-integer, however, this assumption is not valid, and 
charge fluctuations must be treated non-perturbatively. 

To have an insight how change fluctuations can lead 
to non-perturbative behavior, let us study the simplest 
circuit one can envision, the so-called single electron box 
(SEB), where only one lead is attached to a quantum 
dot (see Fig. I19fl . Let us furthermore focus to the limit 
of small tunneling and Eq » T » A. The charging 
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energy of the dot in this case is given again by Eq. J5J, 
and the tunnel coupling to the electrode reads 

3 e 




(a) (b) 

FIG. 19: Sketch of the single electron box. Fig. b. shows 
the top view of the regime where electrons move in the the 
two-dimensional electron gas. Black areas indicate various 
electrodes necessary to shape the electron gas. 

Let us now focus our attention to the regime N g s» 1/2, 
where with a good accuracy the charge on the dot fluc- 
tuates between ndot = and 1. In this regime tunneling 
processes generated by Eq. (|3U[) become correlated, be- 
cause of the constraint that the charge of the dot must 
be either or 1. To keep track of this constraint, we 
can introduce the pseudospin operators, T + = |1)(0|, 
T~ = |0)(1|, and T z = (|1)(1| - |0)(0|)/2, and rewrite 
the tunneling part of the Hamiltonian and -ffdot as 

v = E E{ T+t ^ + h - c -} ' ( 31 ) 

3 e 

fldot = —AS T z (32) 

where AE ~ (N g — 1/2) is simply the energy difference 
between the two charging states of the dot. We can 
now attempt to compute the expectation value of the 
dot charge, (ndot) in the regime —1/2 < N g < 1/2 by 
doing perturbation theory in V to obtain in the limit of 
vanishing level spacing, 

with g = G/Gq the dimensionless conductance of the 
junction. Although a hnite level spacing cuts off the 
logarithmic singularity at N g = ±1/2, Eq. (|33|l clearly 
indicates that perturbation theory breaks down in the 
vicinity of the degeneracy points. 

In fact, following the mapping originally proposed by 
Matveev^4 we shall now show that the Hamiltonian 
above can be mapped to that of the two-channel Kondo 
problem. To perform the mapping, we rewrite the tunnel- 
ing Hamiltonian in a more suggestive way by introducing 
the new fields normalized by the density of states in the 
box and in the lead, gdot and £>i ea d, respectively 

Da = -== V dj, a , C a = — = V *Pe,a , (34) 
VPbox V Plead 



and organize them into a four component spinor 

V^(%). (35) 

In the limit A — > the tunneling amplitudes can be 
simply replaced by their average, tj — > t = (\tj\ 2 )j , and 
we can rewrite the tunneling part of Hamiltonian as 

where the operator just flips the orbital spin r of 
the field tp T<7 , and j± is a dimensionless coupling pro- 
portional to the tunneling, j±_ = 2t A /g> D ox Plead- Thus 
j\ ~ i 2 £?box Plead is directly related to the dimensionless 
conductance g of the tunnel junction 

9 = * 2 j1 ■ (37) 

Eq. Ij36(l is just the Hamiltonian of an anisotropic two- 
channel Kondo modelri^ the orbital spins T and r play- 
ing the role of the spins of the original two-channel Kondo 
model, and the electron spin er providing a silent chan- 
nel index. The presence of this additional channel in- 
dex (electron spin) makes the physics of the two-channel 
Kondo model entirely different from that of the single 
channel Kondo problem, and leads to non-Fermi liquid 
properties^ The low temperature conductance between 
the dot and the lead turns out to scale to g(T — > 0) — > 
e 2 /h with a y/T singularity and the capacitance diverges 
logarithmically as N g -> 1/2 and T -> 0^4 Thus the the- 
ory of Matveev predicts that the sharp steps in (ndot) 
become smeared out as shown in Fig. |3 although the 
slope of the steps diverges at the degeneracy points. 

This is, however, not the full picture. Very recently, 
Karyn le Hur studied the effect of dissipative coupling 
to other leads in the circuit, and showed that if these 
additional leads are resistive enough, then the dissipa- 
tion induced by them leads to a phase transition, where 
the steps are restored^ This transition can be shown 
by bosonization methods to be of Kosterlitz-Thouless 
typei££ In Fig.|20|we show the shape of the step computed 
using numerical renormalization group methods for the 
single channel case (spinless fermions) that clearly shows 
the above phase transition^ 

We emphasize again, that the above mapping holds 
only in the regime A < T, oj, ... < Ec, where the level 
spacing of the box can be neglected, and is only valid 
for the case of a single mode contact, Unfortunately, for 
small semiconductor dots with large enough Eq the ra- 
tio Ec I A is not very largejii*^ and therefore the regime 
where the two-channel Kondo behavior could be observed 
is rather limited. In fact, this intriguing two-channel 
Kondo behavior has never been observed convincingly 
experimentally^ 

The ratio Ec I 'A can be much larger in metallic 
grains, 4,10 However, metallic grains have been connected 
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FIG. 20: The expectation value of ridot — 1/2 for different 
values of a. For moderate values of the bosonic coupling the 
step is smeared out by quantum fluctuations of the charge of 
the dot. The step gets sharper and sharper as a approaches 
a c , and For a > a c a jump appears in (ridot) (from Ref. IHt1> . 



so far only through leads with a large number of conduc- 
tion modes, and the behavior of these systems is rather 
different from what we discussed so far: They can be de- 
scribed by an infinite channel Kondo model;* which pre- 
dicts, e.g., that the junction conductance goes to zero log- 
arithmically at the degeneracy point in the T — > limit, 
9met&\{T) ~ 1/ \n 2 (T/T[<), in clear disagreement with the 
two-channel Kondo predictions. It turns out, that for a, fi- 
nite number N c of conductance modes a new temperature 
scale T* appears*^ Above this scale the conductance de- 
creases as ~ 1/ In 2 (T) , while below this scale it increases 
and approaches the two-channel Kondo value g = e 2 /h. 
Unfortunately, this scale goes to zero exponentially fast 
with the number of modes, T* ~ exp(— C N c ). There- 
fore, one really needs to prepare single mode contacts 
to observe the two-channel Kondo behavior in metallic 
grains, which is a major experimental challenge. 

Recently, another realization of the two-channel Kondo 
behavior has also been proposed in semiconductor de- 
vices, where one has a somewhat better control of the 
Kondo scale and Eq / A than in single electron boxesiSiSi 

So far we discussed only the limit of small level spacing 
A — > 0. The physics of the degeneracy point N g = 1/2 
remains also rather non-trivial in the limit T <C A. In 
this regime the dot behaves as a mixed valence atom, 
and charge fluctuations are hugeJi The description of 
these mixed valence fluctuations is a rather complicated 
problem, and requires non-perturbative methods such as 
the application of numerical renormalization groupS 



VI. CONCLUSIONS 

In the present paper we reviewed some of the interest- 
ing strongly correlated states that can be realized using 
quantum dots. These artificial structures behave in many 
respects like artificial atoms, excepting the difference in 
the energy scales that characterize them. Having a full 



control over these devices opened up the possibility of 
building structures that realize unusual strongly corre- 
lated states like the ones discussed in this paper, that 
are very difficult to observe in atomic physics. This tech- 
nology also enabled one to study out of equilibrium and 
transport properties of such individual 'atoms'. 

However, building quantum dots and quantum wires 
is just the first step towards a new technology, which 
aims to construct devices from real atoms and molecules 
instead of mesoscopic structures. Although a real break- 
through took place in recent years, and molecules have 
been contacted and used to construct single electron 
transistors^ the technology is far from being controlled. 

Building (gated) quantum dot arrays in a controlled 
way has not been solved satisfactorily yet either. Hav- 
ing a handle on such systems could give way to study 
quantum phase transitions in lattices of artificial atoms 
in the laboratory. Experimentalists are also facing the 
challenge of producing hybrid structures at the nanoscale: 
these structures can hopefully be used in future spintron- 
ics applications or quantum computing. 

There is a lot of open questions on the theoretical side 
as well: As we mentioned already, the problem of treating 
strongly correlated systems in out of equilibrium is un- 
solved even for simple toy models, and it is still a dream 
to study molecular transport through correlated out of 
equilibrium atomic clusters by combining ab initio and 
many-body methods. Moreover, many important ques- 
tions like the interplay of ferromagnetism and strong cor- 
relations in ferromagnetic grains have not been studied 
in sufficient detail. 
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APPENDIX A: HARTREE APPROXIMATION 
FOR A QUANTUM DOT 

To derive Eq. J2Jl, let us consider a metallic grain 
within the Hartree approximation. Then the wave func- 
tions ipj must be determined self-consistently by solving 
the following equations: 

V H (r) = V(f) + J dr*U(r- f)g(f) , (Al) 

where V is the confining potential generated by the 
positively charged ions, U is the electron-electron in- 
teraction, and Vh is the Hartree potential. The elec- 
tronic density g(r) in Eq. (|Alfl must be computed as 
q(^) = Tlj a n 3,a l^il 2 ' w ith n j.<r the occupation number 
of the levels. 
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The total energy of the system with N — J2j „ n j,& 
electrons is 

E» = Y, E i n ^~\ J drd78{r)V(r-?) Q {?) , (A2) 

where the second term compensates for overcounting the 
electron-electron interaction. In order to compute the 
energy cost of adding another electron to the grain, we 
should solve selfconsistently Eqs. IjAlfl for £\ g 7ij i(T = 

N + 1, compute E^ +1 and then determine the difference 
between these two energies, E% +1 - Eg. However, one 
can approximately compute this energy by just notic- 
ing that the charge of the extra electron in state N + 1 
must go to the surface of the grain to produce an ap- 
proximately constant potential inside the grain, and is 
screened within a layer of the Fermi wavelength <~ Xp. 
The change in the Hartree potential is simply given as 

SV H = J dr'Uif- r 1 ) 5q(?) . (A3) 

But since the change of the electronic density can be 



very well approximated by a classical surface charge for 
a grain size L ^> Xp, SVh is just the classical potential 
of the charged grain. Consequently, SVh ~ e 2 /C inside 
the metallic grain. Using this simple fact we find that 
adding an extra electron to the grain shifts all Hartree 
energies as Ej — > Ej + e 2 /C and requires an energy 



E+ ~ E f + %n ' ^ A4 ) 



with Ef the Hartree energy of the first unoccupied level, 
and C the classical capacitance of the grain. Defining the 
chemical potential as the Hartree energy of the last occu- 
pied level, /i = e; and defining the quasiparticlc energies 
as ej = Ej — ji we find by extending the above analysis to 
the excited states as well that the energy of the dot is ap- 
proximately described by Eqs. JIJ and Note that for 
an isolated dot this analysis gives N g — —[iC/e 2 , which 
is usually not equal to zero, so adding and removing an 
electron requires different energies, just as in case of an 
atom. 
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